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Abstract 

HST transit observations in the near-UV performed in 2009 made WASP- 12b 
one of the most "mysterious" exoplanets; the system presents an early-ingress, 
which can be explained by the presence of optically thick matter located ahead 
of the planet at a distance of ~4-5 planet radii. This work follows previous 
attempts to explain this asymmetry with an exospheric outflow or a bow shock, 
induced by a planetary magnetic field, and provides a numerical solution of the 
early-ingress, though we did not perform any radiative transfer calculation. We 
performed pure 3D gas dynamic simulations of the plasma interaction between 
WASP-12b and its host star, and describe the flow pattern in the system. In 
particular, we show that the overfilling of the planet's Roche lobe leads to a 
noticeable outflow from the upper atmosphere in the direction of the Li and L 2 
points. Due to the conservation of the angular momentum, the flow to the Li 
point is deflected in the direction of the planet's orbital motion, while the flow 
towards L 2 is deflected in the opposite direction, resulting in a non-axisymmetric 
envelope, surrounding the planet. The supersonic motion of the planet inside the 
stellar wind leads to the formation of a bow shock with a complex shape. The 
existence of the bow shock slows down the outflow through the Li and L 2 points, 
allowing us to consider a long-living flow structure which is in the steady-state. 

Subject headings: hydrodynamics — planet-star interactions — stars: individual 
(WASP-12) 
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Introduction 



WASP-12 is a late F-type m ain sequence star (M « = 1.35 M Q , i?* = 1.6 i? Q , T e g = 
6250 ± 100 if, log g = 4.2 ± 0.2; iFossati et all boiObh wi th a magnitude of V ~ ll m .6. 
The star is at a distance of about 400 pc from the Sun (IFossati et al.l l2010bl ). The star 
hosts a transiting "hot- Jupiter" WASP-12b, wit h a mass of M pl = 1.41 ± 0.1M Jup and 
radius of -R p i = 1.74 ± 0.0 9 -R/i lip (IChan et al.ll201ll ). WA SP- 12b revolves in a rather circular 
orbit (ICampo et al.ll201ll ) with a period of ~ 1.09 days (IHebb et al.ll2009l ) at a distance of 
0.0229 AU (~3 stellar radii) from its host star. 

WASP- 12b was obs erved in 2009 in the near-UV with the Cosmic Origins Spectrograph 
(COS) on board of HST (IFossati et al.ll2010al ). The observed spectral range was 2540-28 10A, 
divided into three bands of 40 A each. Analysis of the COS light curves showed that the 
transit, in two of the three observed spectral regions, is considerably deeper than in the 
visible. Besides, an early ingress was also discovered. The HST data show that in the near- 
UV the transit starts approximately 50 minutes earlier than in the visible, meaning that 
ahead of the planet, at a distance of 4-5 planet radii, there is a region of absorbing matter. 

So far, two explanations for this optically thick matter ahead of the planet have been pro 



posed. The first exp lanation is based on mass transfer from the planet to the star ( iLai et al. 
201CH : iLi et al.ll2010l ). Indeed, the distance of the center of the planet to the Li point is only 
~ 1.85 Rpi, and its upper atmosphere exceeds the Roche lobe. 

The size of the collision dominated atmosphere is limited by the exobase level R ex , where 
the mean free path is equal to the atmosphere's scale height. For WASP-12b's atmosphere, 
which has an expected density at the base of the thermosphere of p p \ = 2.7 x 10~ 14 g/cm 3 
(see Sect. [2]) and is heated by the stellar EUV radiation to a temperature of T p i = 10 4 K, 
the exobase would expand beyond the Roche lobe distan ce. Using standard physics of bi- 



nary stars (see e.g., IPringle fc Waddll985l : ISavonijelll979l ). to estimate how much the upper 



atmosphere overfills the Roche lobe, one can calculate the volume Roche radius R^, ob- 
taining ~ 1.37 i? p i and the degree of overfilling AR/R^ = (i? ex — -RlJ/^Li ~ 0.13. Such 
a large degree of overfilling leads to a powerful outflow of material through the Li and L2 
points. The natural result of this outflow is the formation of an accretion disk surrounding 
the star, as for close bina ries, where the flow struct ure is well studied, both analytically (see, 
e.g., Lubow fc Shul 19751) and numerically (see, e.g.. iBisikalo et al.ll2003c iBisikalo fc Fridman 



2008 



Zhilkin et al 



20121 ) . 



If the accretion disk is axisymmetric it does not produce any effect on the transit shape. 
However, in the re gion where the stream interacts with the disk an observable extended 



shock wave forms (iBoyarchuk et al.l 120021 ; IBisikalo et al.ll2003l ). This hypothesis was used 
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by lLai et al.l ( 120101 ). who supposed that the stream from Li itself and/or the region of the 
stream-disk interaction can cause the early ingress. Although the idea is reasonable, a more 
detailed analysis shows that, with such a large degree of Roche lobe overfilling, the lifetime 
of the planet atmosphere would be very short. In acc ordance with equations used in the 
theory of close binaries (see e.g., iPringle fc Waddll985l ). the mass loss rate is defined as: 



M/M = (AR/Rf 



GM 



Using this equation and assuming a total mass of the WASP-12b's atmosphere of 1.41Mj up , 
we find that the lifetime of the atmosphere is no longer than several years. Accounting for the 
adjustment of the mass loss rate with time, one can then more accurately estimate the lifetime 
of the outflowing atmosphere. Let us assume that the temperature of the th ermosphere is 



const ant, which is a reasonable assumption as WASP- 12 is an inactive star (IFossati et al. 



2010bl ). In this case the mass decrease of a gas giant planet leads to an increase of the Roche 
lobe overfilling, and consequently an increase of the mass loss rate. As a matter of fact, 
for a gas giant planet, equation ([I]) gives an estimate of the maximum planet life time. On 
the other hand, for a planet with a rocky core the mass of its atmosphere is just a small 
fraction of the total mass, therefore one can assume that the size of the Roche lobe will be 
constant, independently of the atmosphere loss rate. Hence, the part of the atmosphere, 
which lies above the Roche lobe, will flow out till the atmosphere itself becomes smaller than 
the Roche lobe, arresting the mass loss. In this way, for rocky planets one can use equation 
(PQ) to obtain a reasonable estimate of the period of time in which the planet is actively losing 
mass. In other words, if the overfilling is high, the lifetime of the atmosphere will be short 
for any kind of pla net, hence a small ch ance to observe it. Therefore the hypothesis proposed 
by iLi et al.l ( 120101 ) and lLai et al.l (120101 ) , most likely does not allow an adequate explanation 
of the processes ongoing in the WASP- 12 system. Nevertheless this probability is not zero; 
it is indeed quite possible that for system s where the dynamic pressure of the wind is small 
the scenario proposed by lLai et al.l (120101 ) is realistic. To some extent our work is a further 
development of this idea. 

The second possible explanation is the presence of an optic ally thick bow shock regio n 
ahead of the planet, which possesses a proper magnetic field (jVidotto et al.l |2010| . 1201 ll ). 
This interpretation is based on the fact that the planet orbits with a supersonic velocity 
inside the st e llar w ind, forming a bow shock at a distance of ~4-5_R p i ahead of the planet. 
Llama et al.l (120111 ) performed radiative transfer calculations based on Vidotto et al.'s model, 
showing that multiple shock configurations would fit the near-UV light curve, hence the need 
of further observations to constrain the model parameters. 



Ionov et al.l ( 120121 ) adopted a model similar to that of IVidotto et al.l (120101 ) . but assumed 
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that the planet does not possess any magnetic field. In this case the bow shock forms imme- 
diately at the level of the upper atmosphere. Although the authors were able to qualitatively 
model the light curve, the shock results to be at a distance significantly different from the 
observed one. The distance between the contact discontinuity and the center of the planet 
can be defined analyti cally considering the dynamic pressure balance from both sides of the 



contact discontinuity (IBaranov fc Krasnobaevlll977l ). For the considered case, this distance 



should not be greater than y 1. 



from its center ( iVerigin et al.ll2003l ). half of that indicated b y the observations (see 



(see 


Lai et al. 


Vidotto et al. 



2010). It is important to note that both bow shock models (jlonov et al.ll2012t IVidotto et al. 
20101 ) do not account for the presence of the outflowing atmosphere. In presence of an over- 
filling Roche lobe the balance equations should take into account the dynamic and thermal 
pressure of the outflowing atmospheric gas. In this way, t he positions of the bo w shock and 
contact discontinuity will be defined by the gas properties. IVidotto et al.l (120101 ) showed that 
the thermal pressure is negligible for a barometric distribution of the atmospheric gas. In 
case instead of an outflowing atmosphere, the gas from the vicinity of the Li point will form 
a stream with a slowly decreasing density, hence increasing the importance of the thermal 
pressure. More importantly, the stream from the Li point accelerates in the gravitational 
field of the star, increasing its radial velocity as y/r. As a consequence, the dynamic pressure 
of the stream will increase (~ r), while the importance of the magnetic field will decrease as 



(5 2 



This means that for an outflowing atmosphere the role of the stream in the 



solution should be extremely relevant. 

The aim of this study is the application of a 3D gas dynamic model to the gaseous 
envelope of WASP-12b, in order to provide an alternative explanation of the observed early 
ingress. This work will be followed by a more comprehensive study describing the sensitivity 
of the results on the assumption and presenting comparisons between the model and the 
available observations. The work is organized as follows: in Sect. 2 we describe our model 
and the assumed input parameters; in Sect. 3 we present the results of the numerical 
simulations; and in Sect. 4 the main results are discussed. 



2. Model description 



The atmospheric parameters of WASP-12b are not well known, therefore we assume 
an isothermal upper atmosphere with a temperature, which is defined by the host star's 
deposition of XUV-energy, and is similar to that obtained for other hot- Jupiters: T p \ = 10 4 K 
( lYelldl2004t iKoskinen et al.ll2012[ ). Note that the assumed temperature of the th ermosphere 
does not correspond to the equilibrium temperature, which is about 2500 K f lHebb et al. 
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20091 ). but to that of the lower thermosphere which is close to th e visual radius and is heated 
by the stellar EUV radiation ffYelle!l2004l ; iKoskinen et aDl2012h . 



We do not include stellar irradiation in the model, because this will affect just the 
boundary temperature value. In our simulation, we attempt to recover the presence of an 
early ingress at the observed distance, by changing the boundary conditions to obtain a 
solution where the stream from the Li point is stopped by the stellar wind at the observed 
distance. On the other hand, the influence of the boundary temperature on the solution is 
rather important, because it defines the density of the atmospheric gas at the point. This 
is a key parameter because it is used to gather the Roche lobe overfilling and the power of 
the stream from the Li point; this consequently defines the solution. We plan to include a 
more accurate description of the irradiation in a future work. 

We estimated the density of the lower thermosphere considering that the optical depth 
along the line of sight is r = n p \ x l p \ x k p \ = 1. Here, the distance l p \ corresponds to the path 
followed by the line of sight through the spherical layer [R p \ — (R p \ + H p \)], where H p \ is the 
scale height. To avoid uncertainties due to the calculation of the opacities fc p i, we assume that 



the at mosphere of WASP-12b is hydrogen dominated, similarly to what IVidal-Madjar et al. 



(120031 ) considere d for HP 209458b, for whi ch the number density at the photometric radius is 



2 x 10 10 cm" 3 flMurrav-Clav et al.ll2009f ). Taking into account that WASP-12b is ~ 2 times 



heavier and ~ 1.6 times larger (radius) than HD 209458b, we estimated the hydrogen number 
density at the photometric radius as n p i w 1.6 x 10 10 cm -3 , corresponding to a density of 
« 2.7 x 10~ 14 g/cm 3 . We defined the obtained dens ity at the photo metric radius, and the 
corresponding exobase at 1.55R p i. For comparison, lLai et al.l (120101 ) define the density at 
the 1 bar level, and their exobase at 1.59R p i, in agreement with our approach. 

We consider a system configuration, comparable to a binary system, consisting of the 
star with M* = 1.35 M and WASP-12b with mass M pl = 1.3 x 1O~ 3 M . Following the 
observations, we assume that the components of this binary system, having an orbital sepa- 
ration A = 4.9 R Q , move in a circular orbit with a period P or b = 26 h . The linear velocity of 
the planet in this system is 230kms _1 . 

The flow is described by a 3D system of gravitational gas dynamic equations closed 
by the equation of state of a perfect neutral monatomic gas. In this model we neglect the 
non-adiabatic processes of radiative heating and cooling. An analysis of our solution shows 
that the density of the envelope is rather high, so everywhere in the envelope the Knudsen 
numbeiQ is < 1. In the important regions of the solution (along the streams from the Li and 



1 The Knudsen number is defined as: Kn = A/ H, where A is the mean free path, equal to l/(n x a) where 
n is the number density and a is the cross-section, a = 1CP 15 cm 2 , and H is a typical scale of the solution 
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L2 points) the Knudsen number is always < 0.1, allowing us to assume that a gas dynamic 
approach is valid for the problem considered here. 

In our model we use the gas dynamic equations, with the assumption that the planetary 
magnetic field is small. Indeed , the synchronization of hot-Jupiters occurs within a few of 
Myr (e.g., IShowman fc Guillotl 120021 . and references therein), while the age of the star is 
about 2Gyr ( Chan et al. 201lh . hence the planet should be tidally locked and therefore its 
magnetic field be weak. There are certainly various further hypothesis for the formation of 
magnetic fields in gas giants (e.g., induced fields), but all of them do imply the presence of 
fields which most likely would n ot protect the expan ding upper atmosphere, moreover in this 
case where the star is inactive ( IFossati et al.ll2010bf ). The presence of the large value of the 
magnetic field also contradicts the observations and stellar wind interaction models which 
reproduce the observed Lyq absorption spe ctra obtained during transits of hot- Jupiters (e.g., 
Ekenback et al.ll2010t lLammer et al.ll201lh . 



To solve the system of gas dynamic equations we use a Roe-Osher TVD scheme of a 
high approximation order with the Einfeldt modification. This numerical method allows one 
to study flows with a significant densit y co ntrast. Further detai ls of the numerical model 
can be found in lBoyarchuk et ajj (120021 ) and lBisikalo et al.l (120031 ). 



Calculations have been performed in a rotating coordinate frame where the force field 
can be described by the Roche potential 



GM, 



pi 



'x 2 + y 2 + z 2 y/ {x — A) 2 + y 2 + z 2 




x - A 



At 



pi 



At* + Mpi 



(2) 



here Q is the angular velocity of the system's rotation. 



The origin of the coordinate system is at the center of the star; the X-axis is directed 
to the planet; the Z-axis coincides with the rotation axis of the system and is perpendicular 
to the orbital plane; the F-axis finalizes the right-handed system. The calculations have 
been performed in a rectangular homogeneous grid. The size of the calculation domain is 
(25 x 20 x 10) R p i and the grid resolution is 464 x 363 x 182 cells. The accepted size of the cell 
of 0.05 R p i allows the investigation of all main flow pattern features in the planet's vicinity. 

The boundary conditions have been set as follows. The upper atmosphere is assumed 
to be isothermal (T p \ = 10 4 K) and in hydrostatic equilibrium, i.e. the gas velocity in the 
upper atmosphere is zero. At the level r = R p \, we set the value of the density at the visual 
radius equal to 2.7 x 10 _14 g/cm 3 . At the initial conditions, the density of the planet upper 



H = n/(dn/dr). 
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atmosphere was defined according to the barometric formula from the adopted boundary 
conditions till the distance where the density becomes less than the wind density. Outside 
this region the co mputational d o main was filled by gas of stellar wind. We treated the stellar 
wind as done by IVidotto et al. The particle number density in the stellar wind has 

been set equal to 5 x 10 6 cm~ 3 (IVidotto et al.ll2010l ). The stellar wind parameters proper of 
WASP-12 are unknown, therefore we used that of the solar wind. Namely, the temperature 
of the wind has been set equal t o that of the Sun at the corresponding distance from the 
star T = 10 6 K (jWithbrod Il988l ) . The velocity of the wind is assumed to be 100 km s -1 , 
corre sponding to the velocity of the solar wind at the distance of WASP-12b to the host 
star (jWithbrod 119881 ) . As the wind acceleration mechanism is still an open question, in the 
model we imply the existence of a wind acceleration mechanism which works similarly to 
light pressure. To take this into account we modify the expression for the Roche potential, 
given in Eq. (2), as follows 



CM. 



GM, 



P i 



yjx 2 + y 2 + z 2 sj(x - A) 2 + y 2 




A- 



M, 



v i 



P i 



+ y z 



where the T coefficient is used to account for the stellar wind acceleration and it is set to 
zero in the regions filled with wind matter, and one in the rest of the calculation domain. 
This allows us to avoid the non-physical wind deceleration due to stellar gravity. The proper 
stellar wind velocity is subsonic with a Mach number M = 0.85. However, taking into 
account the supersonic orbital motion of the planet (M = 1.97) the total velocity of the 
planet with respect to the stellar wind is indeed supersonic with the rather large Mach 
number of M — 2.14. 



3. The structure of the planet gaseous envelope 

The general morphology of the flow is shown in Fig. [TJ indicating the density distribution 
and velocity vectors in the planet's envelope. The planet is depicted by the filled circle and 
it moves counterclockwise. The solid yellow lines denote the flow lines starting from the gas 
of the stellar wind. Besides, the Roche lobe, depicted by the white solid line, and Lagrangian 
points Li and L2 are also shown. 

Considering the density distribution shown in Fig.[TJ we notice that the planet's envelope 
has a non-spherical complex shape. In addition to the upper atmosphere itself, one can see 
flows towards the Li and L 2 points. According to the conservation of angular momentum 
these streams are deflected in the direction of the planet's motion and against it, respectively. 
The supersonic motion of the planet and the presence of the stellar wind lead to the formation 
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Fig. 1. — Density distribution and velocity vectors in the envelope of WASP- 12b. The planet 
is depicted by the filled circle and it moves counterclockwise. The solid yellow lines denote 
the flow lines starting in the gas of the stellar wind. The white solid lines denote the Roche 
equipotentials passing through the Lagrangian points Li and L 2 . 
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of a bow shock oriented according to the total vector of the wind material velocity with 
respect to the planet. The location of the shock and contact discontinuity can be easily 
determined from variations in the flow lines. It turns out that the head-on collision point is 
located on the gas stream from the L x point. The distance between the planet center and the 
head-on collision point, as projected onto the stellar limb, is ~4.5 planet radii, in agreement 



with the estimation by lLai et al.l ( 120101 ). who derived the position of the early-ingress at 



~3.2 planet radii, from the planet surface. 

Both shock and contact discontinuity have complex shapes. The asymmetric shape of 
the planet envelope, the upper atmosphere itself and two flows/prominences from Li and L2, 
lead to the formation of a distinguishable double-peaked shock. The head-on collision point 
is at the peak of the "prominence" in the direction of Li. However, when moving closer to 
the planet, we see that the shock bends due to the presence of the planetary atmosphere, 
and shifts farther from the planet, leading to the formation of the second "hump" in the 
shock. It is noticeable that the matter's flow, moving to the planet from the head-on collision 
point, undergoes strong disturbance when it gets into the cavity between the two humps in 
the shock waves. In particular, vortices form in this cavity. This effect smooths the contact 
discontinuity, mixing the matter of the stellar wind and upper atmosphere. 

The flows towards the Li and L2 points extend far away from the planet's Roche lobe. 
Formally, the planet loses its atmosphere and, according to (JTJ), the lifetime of the atmosphere 
is not longer than a few years. However, in our simulation the envelope is in a quasi-stationary 
regime, i.e. it is almost enclosed. The dynamic pressure of the stellar wind and flow, caused 
by the orbital motion of the planet, breaks the propagation of the streams from the Li 
and L 2 points, limiting the gaseous envelope of the planet by the bow shock and contact 
discontinuity. 

Let us consider the problem of the system closure in detail. From the physical point of 
view, the material that has left the planet's exosphere through the vicinity of the I_q point 
should follow a certain trajectory determined by the balance of three forces: the force of 
inertia, gravity forces, and the gradient of gas pressure. We can neglect the gravity of the 
planet out of its Roche lobe, hence, the gas motion can be considered as that falling in the 
gravity field of the star. The stream is deflected by the Coriolis force. However, at moderate 
distances from the Li point we can estimate the radial velocity of matter using the difference 
between the values of the potential energy at the point of our interest and the Li point 

v 2 « $(Li) - $(r), (3) 



where $(r) is the Roche potential at the point with the radius- vector r, and Li is the radius- 
vector of the inner Lagrangian point. With the motion in the gravity field of the star the 
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velocity of the stream grows. However, at the same time the wind density also grows in 
accordance with: 

Pro = Pm (J^) ■ (4) 

It results that eventually the dynamic pressures of the stream p s v 2 and wind pwV 2 ^ become 
equal and the radial motion of the stream ends. We only need to estimate at what distance 
this is happening and if this distance is in agreement with the observations. 



According to ballistic analysis ( jLubow fc Shu! Il975[ ) the stream, moving in a binary 



system, is deflected from the line connecting the centers of mass of the star and planet at an 
angle of ~ 20°. The dynamic influence of the wind should increase this angle, and, as follows 
from the results of the numerical simulations (see Fig. [TJ), the angle grows up to ~ 40°. For 
the location of the shock wave to correspond to the observed one (Y w w 5 R p \ from the center) 
the radial motion of the stream should stop at a distance of X w = —5 R p itg40 o ps —6 R p \ 
from the point. We can therefore write the following equation, s etting the equi l ibrium 



between the dynamic pressures of the stream and wind (as done by IVidotto et al.l ( 120101 ) 
but for an atmosphere with a magnetic field), and and (@J 



A^ 2 



p il ($(L 1 )-$(r)) = p W0 ^-j vi. (5) 

Solving this equation we find analytically the density at the Li point whose value allows us 
to stop the shock wave at the observed distance from the planet. Substituting the known 
value of r (from the observations) into (JSj), we obtain wlx 10 _17 g/cm 3 corresponding 
to a number density of ~ 6 x 10 6 cm~ 3 . The obtained value is ~ 18 times less than 
that at the Li point, given by the barometric formula^! where p v \ = 2.7 x 10~ 14 g/cm 3 . In 
the gas dynamic simulation the stream stopped at the observed distance with the density 
value p p i ~ 1.35 x 10~ 15 , which is ~ 20 times less than the boundary density, estimated 
using the planet's photometric radius. In our analysis we omitted the effects of the stream 
gas heating due to the bow shock's radiation, we can state that the analytical estimates are 
in good agreement with the results of the numerical simulations. 

The numerical simulations presented in this work allow us to simulate the formation of 
an enclosed and steady-state gaseous envelope surrounding "hot Jupiters" which exceed the 
planet's Roche lobe, such as WASP-12b. We should note one more i mportant property o f 



the described model of the planet's envelope. Observations indicate ( jFossati et al.ll2010al ) 



that the eclipse in the spectral bands where the early ingress has been observed is two times 



2 Whcn determining this value, we assumed that the density is constant over the whole surface of the 
Roche lobe and is equal to the value found for the shortest distance from the planet to the lobe. 
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deeper than the eclipse of the planet itself (3.2% versus 1.7%). The possible explanation for 
this effect can also be the formation of a bow shock. Indeed, the motion of the planet inside 
the stellar wind is significantly supersonic and the formed bow shock is much hotter than the 
planet upper atmosphere. The heating of the gaseous envelope including the stream from 
the Li point leads to the excitat i on an d bro adening of additi onal lines in the UV bands. 
According to IVidal-Madjar et all (120041 ) and iBen-Jaffell (120071 ) the depth of the transit at 
some spectral bands strongly depends on its width and presence of strong spectral lines. The 
larger the equivalent width of the lines in a certain band, the deeper the transit in this band. 
Therefore, the additional heating of the planet envelope can lead to the observed effect of a 
stronger absorption in the NUV bands. 



4. Conclusions 

Using 3D numerical simulations, we investigated the flow pattern in the vicinity of the 
exoplanet WASP-12b, which is overfilling its Roche lobe. Taking into account that the planet 
is tidally locked we assumed no (or negligible) surface magnetic field and consider the pure 
gas dynamic solution. In this work, we study the flow structure of the forming gaseous 
envelope surrounding the planet, without including any radiative transfer calculation. The 
obtained solution depends strongly on the stellar wind parameters which we assumed close 
to solar. 

Our results indicate that : 

• because the planet's upper atmosphere overfills its Roche lobe, WASP-12b's envelope 
has a complex non-spherical shape. In addition to the central part of the spherical 
upper atmosphere two "prominences", directed to the Li and L 2 Lagrangian points, 
develop in the system. These streams leave the planet and are deflected in the direction 
of the orbital motion and against it. Under the action of the dynamic pressure of the 
stellar wind these flows slow down and, then stop at distances of ~ 6 and ~ 4 planetary 
radii from its center, forming a stationary envelope. 

• The planet and its envelope move in the gas of the stellar wind with a supersonic 
velocity resulting in a Mach number M = 2.14. Thus, the dynamic pressure of the 
stellar wind not only works towards the formation of the stationary envelope, but 
leads also to the formation of a bow shock and contact discontinuity, which delimit the 
envelope. The wave has a complex double-peaked shape. 

• The head-on collision point found in the calculations is at the peak of the stream from 
the Li point. The distance from the planet to the head-on collision point, as projected 
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to the stellar limb, is of ~ 4.5 planetary radii, allowing us to explain the observed 
extent of the early ingress. 

• The heating of the planet envelope by the bow shock allows one to explain the fact 
that the eclipse, in the spectral bands where the early ingress has been observed, is 
two times deeper than the eclipse of the planet itself. 

Summarizing our results, we found that by using the full gas dynamic simulations of 
the stellar wind interaction between WASP- 12b and its host star, we have been able to 
obtain an alternative and self-consistent flow pattern allowing the explanation of the existing 
observational data. 
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